May 2006 



Quantum Fields on Star Graphs 



B. Bellazzini al and M. Mintchev 62 

a INFN and Scuola Normale Superiore, Piazza dei Cavalieri 7, 56127 Pisa, Italy 

b INFN and Dipartimento di Fisica, Universitd di Pisa, Largo Pontecorvo 3, 56127 

Pisa, Italy 



Abstract 

We construct canonical quantum fields which propagate on a star graph 
modeling a quantum wire. The construction uses a deformation of the algebra 
of canonical commutation relations, encoding the interaction in the vertex of 
the graph. We discuss in this framework the Casimir effect and derive the 
correction to the Stefan-Boltzmann law induced by the vertex interaction. We 
also generalize the algebraic setting in order to cover systems with integrable 
bulk interactions and solve the quantum non-linear Schrodinger model on a 
star graph. 
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1 Introduction 

The physics of thin three-dimensional graph-like structures , two of whose dimen- 
sions being of the order of a few nanometers, attracts recently much attention. Due 
to the impressive progress in nanotechnology such devices, called quantum wires, 
can be constructed and tested nowadays in laboratory At the theoretical side, 
the quantum behavior of these essentially one-dimensional systems poses interesting 
mathematical and physical problems, which have been addressed in the past decade 
by various authors In the present paper we pursue further the study of quan- 

tum wires, developing field theory on such structures. For this purpose we adapt 
the algebraic approach developed in ^0]- for dealing with systems with bound- 
aries and point-like defects. This approach translates the boundary value problem 
at hand in algebraic terms and offers an efficient framework for the construction of 
quantum fields. 




Figure 1: A star graph G with n edges. 

We focus in this paper on quantum wires of the form of a star graph G shown 
in Fig. 1, where V stands for the vertex of the graph and (with i = l,...,n) 
denote its edges. Being the fundamental building blocks for general graphs, the star 
graphs play a distinguished role [HUE] in the subject. The first question, we address 
in section 2, is the explicit construction of a canonical scalar quantum field, freely 
propagating in the bulk B = G \ V of the graph G. The vertex is considered in 
our approach as a defect (boundary in the case n — 1) interacting with the field 
and characterized by reflection {n > 1) and transmission (n > 2) coefficients. We 
consider the case without dissipation and treat both relativistic and non-relativistic 
dispersion relations. In section 3 we derive the correlation functions with respect to 
two different cyclic states - the Fock vacuum and a Gibbs state at finite (inverse) tem- 
perature (3. Some basic physical observables are analyzed in section 4. In particular, 
we compute there the Casimir energy and the correction to the Stefan-Boltzmann 
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law generated by the interaction localized in the vertex of the graph. Inspired by 
factorized scattering in 1+1 dimensional integrable systems, we introduce in section 
5 a two-body interaction in the bulk B. The system treated here is the non-linear 
Schrodinger model (non-relativistic ip A theory) on a star graph. We establish the 
exact operator solution of this model. Section 6 contains our conclusions and some 
ideas for future developments. 



2 Quantum fields on a star graph 

Consider a star graph G. Each point P in the bulk B of G can be parametrized 
by a pair (i,x) where i = 1, ...,n indicates the edge and x > the distance of P 
from the vertex V along the edge. We will see below that this choice of coordinates 
ensures a uniform treatment of all edges as far as the reflection and transmission 
by the vertex V is concerned. The embedding of G and the relative position of its 
edges Ei in the ambient space are irrelevant in what follows. 

2.1 Non-relativistic dispersion 

We start with the non-relativistic complex scalar field tp(t,x,i) satisfying 

(id t + ^-d^j^{t,x) = 0, teR, x>0, i = l,...,n, (2.1) 

with standard initial conditions given by the equal-time canonical commutation 
relations 

[V> fl (0, xt) , i> i2 (0, x 2 )\ = {r h (0, xt) , r 12 (0, z 2 )] = , (2.2) 

[i> il (0,x 1 ) , <T 2 (0,x 2 )] = 5tl 8{x 1 - x 2 ) , (2.3) 

where * denotes Hermitian conjugation. Since eqs. (|2.1H2.3|) hold only in the bulk 
B, they do not fix uniquely the solution. For this purpose one needs some boundary 
conditions in the vertex V. Setting 

if)i(t,0) = lim^(t,x), (d x i(>i)(t,0) = ]im(d x il>i)(t,x), (2.4) 

x>0 x>0 

we require jS] that 3 

A?ip j (t,0)+Bi(d x <tp j )(t,0) = 0, Vt G 1, z=l,...,n, (2.5) 

3 Summation over repeated upper and lower indices is understood, unless otherwise stated. 
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A and B being two complex n x n matrices. It has been shown in that (J2.5)) 
uniquely defines a self-adjoint extension of the operator —d 2 x on G, provided that 
the composite matrix (A, B) has a rank n and 

AB*-BA* = 0. (2.6) 

Moreover, this extension admits a non-trivial scattering matrix, reading in momen- 
tum space 

S(k) = -(A + ikB)-\A-ikB), (2.7) 
which, in view of (12. can be written also in the form 

S(k) = -{A* - \kB*){AA* + k 2 BB*)-\A - ikB) . (2.8) 

The entries of S have a simple physical interpretation: the diagonal element Sl(k) is 
the reflection amplitude on the edge E iy whereas S J { (k) with i ^ j is the transmission 
amplitude from Ei to Ej. From (j2.7J) and ()2.8|) one can easily deduce that S is unitary 

S(k)* = S(k)~ 1 , (2.9) 

and satisfies Hermitian analyticity 

S{k)* = S{-k). (2.10) 

Combining ()2.9|) and (j2.10|) one gets 

S(k)S(-k)=I n , (2.11) 

where I n is the n x n identity matrix. The property (j2.11j) is essential in what 
follows. 

Turning back to the quantum field ip, we introduce the wave functions [3] 

xftx; k) = ^ kx 5{ + Sl{-k)e~' lkx , k > , (2.12) 
which are orthogonal (see the appendix) 

dx X *!(x;k)xi(x;p) = 5{27T5(k-p). (2.13) 



If in addition 

-e ife ^(fc) = 0, x>0, (2.14) 
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then x form a complete set, i.e. 



dk 

^Xl(x;k)xi(y,k) = 6l6(x-y). (2.15) 



At this point we define 



A(t,x) = / -e-^xKx; k) aj (k) , (2.16) 
Jo 271 

— e ^a*i(k)x*;(x;k), (2.17) 

where u(k) is the dispersion relation 

k 2 

u(k) = — , m>0, (2.18) 
2m 

and {di(k), a* l (k) : k > 0} generate an associative algebra ^4+ with identity element 
1 and satisfy the canonical commutation relations 

a h (ki) a i2 {k 2 ) - a i2 (k 2 ) a^(ki) = 0, (2.19) 
a* h (k{) a" 2 (Afc) - a" 2 (Jfe 2 ) a" 1 (A*) = , (2.20) 
a h (h) a" 2 (fc 2 ) - a« 2 (& 2 ) (A*) = 2tt5(A; 1 - fc 2 ) $ 1 . (2.21) 

The equation of motion ()2.1|) and the boundary condition (|2.5jl follow directly from 
the representation (J2.16JI . The commutation relations ()2.2l2.3j) are a consequence 
of (j2.19H2.2T|) and the condition (|2.14|) . If the latter is violated, in addition to 
{a,i(k), a* l (k) : k > 0} the system admits other degrees of freedom, which must be 
taken into account for reproducing (|2.Hj) . These new degrees correspond to vertex 
bound states, which have been described in detail for n = 2 in ^3]- In order to 
simplify the discussion, we assume in the present paper that ()2.14j) holds. 

In the representation ()2.16I2.17|) the boundary condition ()2.5|) is captured by 
the wave functions x an d X* ■ Remarkably enough, this information can be shifted 
to a,i{k) and a n (k), extending A + to any k e R and considering the algebra A 
generated by {a^k), a*\k) : k G M} which satisfy (12.1 912. 2()|1 . 

a h (h) a" 2 (k 2 ) - a" 2 (k 2 ) a h (h) = 2n5(h - k 2 ) tfg 1 + S% (k 1 )2n5(k 1 + k 2 )l, (2.22) 

and the constraints 

cn(k) = Si(k) aj (-k) , a H (k) = a^i-^SU-k) . (2.23) 
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One can easily verify in fact that in terms of A the fields if) and ifi* take the simpler 
form 

/ + °° rlL 
^ ai {k)^ {k)t+ikx , (2-24) 

r+°° (ih 

if,«(t,x)= —a*\ky^ k)t -' kx , (2.25) 



— oo 



where now the integration over k runs over the whole line. The constraints ()2.23|) 
implement in ()2.24I2.25|) the boundary conditions ()2.5|) and are consistent because 
of (|2.1ip . Hermitian analyticity (|2.10j) implies that the mapping 

I : a H (k) h-> di(k) , I : a^k) i-> a*%k) , (2.26) 

generates an involution in A, which is consistent with the second term in the right 
hand side of ()2.22|) . Because of f)2. 14j) . the initial conditions (J2.2J2.3)) are not in- 
fluenced by this term. It contributes however to the time evolution and codifies 
the interaction in the vertex V, playing the role of a defect. Indeed, the algebra A 
defined by ()2. 1912.2012.221 12.23|) is a special case of the boundary and the reflection- 
transmission algebras introduced in ^U] and [TT] IT2*] for dealing respectively with 
boundaries and defects in quantum field theory. With our present choice of co- 
ordinates, the direction of both reflected and transmitted waves coincides with the 
orientation of the edges. For this reason both reflection and transmission amplitudes 
in ()2.22|) multiply the 5(ki + /^-contribution. 4 

The Hamiltonian H of our system is the self-adjoint extension of —0%, defined by 
the pair (A, B). The work with this somehow implicit form of H is greatly simplified 
in the algebraic setting, where H takes the familiar quadratic form 

1 f+OO 77 /-+00 77 

H = - —ujikX'i^diik) = / — uj{k)a H {k) ai {k) , (2.27) 



2 J-oQ 27T J 2tt 

which is very convenient to deal with. One easily verifies for instance 

^i(t,a;) = e^ i (0,a;)e- ltH , (2.28) 

which confirms that ()2.27j) generates indeed the time evolution of if). It is also worth 
stressing that the energy is conserved, showing that the vertex V behaves like a 
defect without dissipation. 

We are interested in the fields if> and if)* in the range x > 0. Nevertheless, 
eqs. (j2.24l2.25J) keep a well-defined meaning for any x e R. The resulting fields 

4 This is not the case for the coordinates adopted for n = 2 in ^21> where the transmission 
coefficients multiply S(ki — k^] because E\ = R+ but E<x = M_. 
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satisfy the equation of motion ([2.1)1 . but do not obey the equal-time commutator 
([2.3)1 on the whole line. The constraints ([2.23)1 imply 

&(t,x)= f dySi(y)iP j (t,y-x), (2.29) 

J — oo 

which relates the values of ip on R + with those on M_. 

Summarizing, we translated so far the construction [3| of the self-adjoint exten- 
sion of the operator —d%, relative to the boundary condition ([2.5)1 . in terms of the 
algebra A. The advantage of the algebraic approach is that: 

(i) it admits a straightforward extension to more general scattering matrices then 
those defined by ([2.7)1 : 

(ii) it provides a framework for introducing and studying some integrable inter- 
actions in the bulk of G. 

Let us concentrate now on point (i), dedicating to (ii) section 5 below. The key 
observation is that the algebra A is actually well-defined for any S'-matrix satisfying 
(12. 912. 10)1 . which allows to consider much more general boundary conditions then 
(12.5)) . This fact is evident already for n = 1. In this case the general S- matrix 
satisfying ([2.912. 10)) can be written in the form 

S(k) = - , s(-k) = -s(k) , s(k) G R . (2.30) 

s{k) + i 

Under mild technical assumptions on s, one can define the pseudo-differential oper- 
ator 



s(k)f(k)e- ik % (2.31) 



f being the Fourier transform of /. The function s is called the symbol of the 
operator D s . Setting e p (x) = e~ ipx , one immediately verifies that 

[D s e p ] (x) = s(p)e p (x) , [D s e_ p ] (x) = s(-p)e~ p (x) = -s(p)e^ p (x) . (2.32) 

Using these simple identities and the representation ([2.16)) with S given by ()2.3())) . 
one obtains that ip = ip% satisfies the general boundary condition 

[(£>, + i) V] (*, 0) = . (2.33) 

The above argument extends to star graphs with n > 1 edges as well. One con- 
cludes therefore that besides ()2.5)) . the algebra A covers also more general boundary 
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conditions, formulated in terms of pseudo-differential operators. Computing the cor- 
relation functions of ip and ip* in the next section, we apply the general framework 
offered by A without referring necessarily to ()2.5|) . 

It is well-known that the non-relativistic field ip can be quantized with Fermi 
statistics as well, replacing the commutators (|2.2I2.3|) with anticommutators. Obvi- 
ously, for this purpose one has to perform the same modification in (|2.19H2.2'2j) . 



2.2 Relativistic dispersion 

The above setting has a straightforward extension to quantum fields with relativistic 
dispersion relation. Let us consider for example the Hermitian scalar field tp(t,x,i) 
defined on the graph G by 

(5f - d 2 x + m 2 ) ip{t,x,i) = 0, teR, x>0, i = l,...,n, (2.34) 

and the equal-time canonical commutation relations 

[<p(0, xx,i t ) , tp(0, x 2 , i 2 )] = , (2.35) 

[(3tf>)(0, x 1 ,i 1 ) , <p(Q, x 2 , i2)\ = -Mil S(x, - x 2 ) . (2.36) 

Since ip is Hermitian, we impose on it the counterpart of the vertex boundary con- 
dition ()2.5|) for real A and B. In this case one infers from (J2.8|) that S is symmetric, 
i. e. 

S\k) = S{k) . (2.37) 

Like in the non-relativistic case, we introduce at this stage the algebra A with S 
obeying (I19I2U()I114I2\37I) . Then 

j- [a* l {k)e iu)ik)t - [kx + ai{k)e- iu){k)t+ikx ] , (2.38) 

-oo 27Ta/2u;(/c) 

with 

oj(k) = Vk 2 + m? , (2.39) 

is the solution of (|2.34H2.3^j) we are looking for. The Hamiltonian is obtained from 
(|2.27|l . inserting the relativistic dispersion relation ()2.39|) . 

Summarizing, we introduced canonical quantum fields on a star graph G. The 
new feature, with respect to the conventional free fields on M s , is the presence of 
a non-trivial one-particle S'-matrix describing the interaction at the vertex V. The 
latter represents physically a point-like defect. Quantum field theory with such 
defects has been investigated |IS]-[2ni i n various frameworks. The results of this 
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section confirm once more the conclusion of [21], indicating the algebra A as a 
universal tool for quantization with defects. Our goal in what follows will be to 
investigate the basic features of the quantum fields constructed so far on G. Before 
doing that however, it is instructive to illustrate at this stage the vertex boundary 
conditions ()2.5|) with some examples. 



2.3 Examples 

The simplest choice for A and B in ()2.5|) is 



/ 1 -1 
1 -1 



.4 




\0 

which give rise to the star graph generalization 

Vi(*,0)=Va(t,0) = ...= 



\ 




1 -1 
1 



B = 



^„(t,0) =0 



of the familiar Dirichlet boundary condition. In this case S(k) = I n . 
A slightly more complicated example is 



/ 1 -1 
1 -1 



A 








\ 



1 -1 

o -v J 



B 



( 





1 1 1 



\ 





1 1 



which leads to 

^i(t,0)=^ 2 (t,0) 



i=l 



;2.40) 



(2.4i; 



;2.42) 



VnM), ^T l {d x >i/> i )(t,Q)=rfil> n {t,0), (2.43) 



generalizing the Robin (mixed) boundary condition. The associated S-matrix is 
non-trivial, 



S(k) 



( (2-n)k-irj 2k 2k ■ ■ ■ 

2k (2-n)k-irj 2k ■■■ 



nk + irj 



2k 
2k 



\ 



2k 



2k 



2k 



(2 — n)k — ir) j 



(2.44) 
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In particular, 77 = provides a generalization of the Neumann condition for n > 2. 

Notice that the boundary conditions (|2.41I2.43|) are symmetric under edge per- 
mutations. This is clearly not the case in general. A simple asymmetric example is 
defined by 




-1 


1 








-1 


1 












.4 



The associated ^-matrix reads 



A > 



B 




(2.45) 



S(k) 



3(k + i\) 




(2.46) 



which shows that the boundary condition determined by ()2.45|) is not invariant under 
the permutations 1 2 and 1 «-> 3. 

Being associated with critical points in the context of statistical mechanics, the 
scale invariant vertex boundary conditions are of special interest. The subset of 
such conditions of the form (j2.5J) has been fully described in |3J. As expected, the 
corresponding S'-matrices are k- independent. An example of this type, obtained by 
setting 7] = in ()2.44|) . is 



S 



( (2-n) 2 2 • 
2 (2-n) 2 • 



n 



\ 



\ 



2 • 



(2.47) 



■ (2-n) J 



Let us observe finally that given a pair of matrices (A, B) fixing the boundary 
condition (|2.5jl . the pair (A,B) = (—B,A) defines an admissible condition as well. 
By means of (|2.8|) . the relative scattering matrix satisfies 



S(k) = S(k~ 



(2.48) 



This remarkable property e 3 3 shows that the mapping (A, B) i— > (-B, A) is a sort 
of duality transformation, relating high and low energies. Being fc-independent, the 
scale invariant S'-matrices are obviously dual invariant. One should keep in mind 
that duality does not preserve in general the completeness condition (|2.14jl . In fact, 
it follows from ()2.48|) that resonant states are dual images of bound states and vice 
versa. 
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3 Correlation functions 

In order to extract concrete physical information from the algebra A, one should 
study its representations. We focus here on two of them: the Fock representation 
J 7 (A) and the Gibbs representation Gp{A) at inverse temperature (3. Since all tech- 
nical details about the construction and the properties of J 7 (A) and Gp(A) can be 
found in and respectively, we collect below only the basic formulae. Em- 
ploying the commutation relations (|2. 1912. 2012. 22j) . one can reduce the computation 
of a generic correlation function to correlators of the form 

m n 

fc=i i=i 

which can be evaluated in turn by iteration via 

m n m m n 

k=l 1=1 k=l k=2 1=1 

l^k 

(3.2) 

For this purpose one needs only the two-point functions. In J 7 (A) on has 

(ai( P )a*i(q)) = 2vr [5{ 5(p - q) + S*(p) 5(p + q)] , (a* i (p)a j (q)) = . (3.3) 
In Q 13(A) one finds ^3] instead 

(^(PK^)^ = l ± e „l Mp) ^ } 2vr [5iS(p -q) + Sj(p)S(p + q)} , (3.4) 

(a*\p)aM))p = l ± e ^ Mp) ^ } 2^ [8)8{p -q) + S}(-p)5(p + q)] , (3.5) 

where ± stands for Fermi/Bose statistics and fi G R is the chemical potential asso- 
ciated with the number operator 

N=-J ^ -a»(k) ai (k). (3.6) 

A characteristic property of any thermal representation is the Kubo-Martin-Schwinger 
(KMS) condition. In our context it states that 

([a a (n(p)] a* J (q)} p = (a* J (q) [a s+i/3 ai(p)})p , (3.7) 

a s being the automorphism 

a s dT{k) = e^a^ikje-^ , a s a^k) = e isK ai {k)e- isK , K = H - fiN . 

(3.8) 
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One easily verifies that (|3.4I3.5|) satisfy (|3.7|) . 

It is straightforward now to derive the two-point correlators of the fields ip, ip* 
and if introduced in the previous section. In the Fock representation J-"(A) one gets 

(1> il (t 1 ,x 1 )i> ia (t 2 ,x 2 )) = (r n {ti,x 1 )r i2 {t2,x 2 )) = 0, (3.9) 

/ + °° AU 

_ e -Mk)t 12 ^ 12J j + e i^ 12 ^ (fc) ] ; (3 . 10 ) 

with 

x 12 = x 1 - x 2 , x± 2 = x 1 + x 2 (3-11) 
and the dispersion relation (j2.18j) . In the relativistic case (I2.39j) one obtains 

/+°° au 
i ^e-^-[e i ^Mj + e ife ^^(A ; )] . (3.12) 

We see that the fields in the different edges are not independent, but interact through 
the scattering matrix S taking into account the boundary conditions in the vertex 
of the graph. 

In the Gibbs representation Gp{A) one finds 5 

(^ il {t 1 ,x 1 )ij; i2 {h,x 2 )) p = (^{tiiX^ihiX^p = 0, (3.13) 

^ il {ti,x l w i2 {^x 2 ))p = 

/+°° AU 1 
— m m , e- [uj{k)t12 \^ Xl2 ST + e ikx ™S?Ak)\ , (3.14) 

(4>* 11 (h, xi)ij i2 (t 2 , x 2 )) i3 = 

r+°° AU p -/3{u(k)-n] 

((p(t 1 ,X 1 ,i 1 )(p(t2,X2,i 2 )}{3 = 

/ + °° dk p-/3[^(fe)-At]+iw(fc)ti2 _i_ p -iLu(k)ti2 
„ ^J(k) ^M^S ^ + e^(*)] . (3.16) 

With the above background we are ready to compute some quantum field theory 
observables like the particle, current and energy densities on a star graph. 



5 In the bosonic case we assume that /i < 0. 
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4 Observables 

The basic local observables associated with the field ip are the particle density 

g{t,x,i)= [r%](t,x), (4.1) 

the current density 

j{t,x,i) = ^ [rXd^i) - [d^A] {t,x) , (4.2) 
and the energy density 

0(t, x, i) = ~ ty™ (d 2 x A) + {dlr 1 ) A] (t, x) , (4.3) 
without summation over % in the right hand side. Conservation and Kirchoff 's laws 

n 

d t g{t,x,i) = d x j(t,x,i), ^2j(t,0,i) = 0, (4-4) 

t=i 

follow from (|2.1J) and (|2.5jl respectively. 

It is straightforward to compute the expectation values of (j4.lH4.3j) in the Gibbs 
representation Gp(A) defined in the previous section. Applying the conventional 
point-splitting procedure, one gets from (|3.15|) 



(g(t,x,i))p = lim (ip* ll (t 1: Xi)ipi 2 (h, x 2 )}/3 

x^ — >X2 = a; 

"+°° dk e -/?Kk)-ri ^ k 2 



Analogously, (j(t,x,i))p = and 

/■ +0 ° r/fr p-|8M*)-f] 



Egs. 1)4.514.6(1 provide convenient integral representations of the particle and energy 
densities in the Gibbs state. Since the integrands involve only diagonal elements 
of the S'-matrix, the expectation values (|4. 514.6(1 are real because of Hermitian an- 
alyticity ()2.10j) . In agreement with the invariance of the Gibbs state under time 
translations, (j4.5l4.6j) are t-independent and vanish in the limit of zero temperature 
((3 — * oo), which shows that these densities have a purely thermal origin. 
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Let us concentrate on the energy density (|4.6|) . It is instructive to separate the 
universal contribution 



from the vertex dependent part 

Unitarity ()2.9[) of S implies the estimate 

|£±0M,/3)| <e±{/3), (4.9) 

leading in turn to positivity 

(9(t, x, i)) p = e ± {(5) + S ± (x, i, (3) > (4.10) 
of the energy density in the Gibbs state. The integration over k in ()4.7|) gives 

3 

e ±(0) = T^T^Lis ( T ^) f ^ * , (4.11) 



where Li s is the polylogarithm function. Eq. (|4.11j) is the analog of the Stefan- 
Boltzmann (S-B) law for the non-relativistic dispersion relation (['2.180 . With /i = 
it simplifies further 



C being Riemann's zeta function. 

According to the previous discussion, the vertex dependent contribution (J4.8[) 
represents a correction to the S-B law generated by the defect at x = 0. As ex- 
pected, this term is both x and /3-dependent and involves the S-matrix. The explicit 
integration in k is usually a hard problem, but the integral representation (j4.8j) is 
very convenient for the numerical evaluation of £±(x, i, (3). 

Let us consider now the case with relativistic dispersion (j2.39j) . focusing on the 
scalar field (J2.34H2.3?)|) . The energy density in this case reads 



6{t, x,i) = -\ (%?) 2 - - [<p{%<p) + {ftp)<p] + mV 2 (t, x, i) . (4.13) 
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One can compute the expectation value (9(t, x, i))p via point-splitting from the two- 
point function (J3.16j) . Differently from the non-relativistic case, it turns out that 
(6(t,x,i))j3 diverges. This is not surprising because the same problem is encoun- 
tered already in the standard Casimir effect. Like in that case, one can solve it by 
subtracting from (9(t, x, the expectation value (9{t, x)) 1 ™ of the energy density 
8(t, x) on the whole line R and at zero temperature. The result is 

(9(t, x, i)) p - (6{t, x)) 1 ™ = e s _ B (P) + £ c (x, i) + £(x, i, (3) , (4.14) 

where 

,+ °° dk „ N e-PWQ-A 



£ S-B 



W = ^^) 1 _ e -^( fc )-,] > "(*) = v^T^, (4.15) 



is the Stefan-Boltzmann contribution, 

1 f + °° dk 

€ a (x,i) = - J ^u(k)SKk)e 2 ^ (4.16) 

is the Casimir energy density associated with the vertex interaction at zero temper- 
ature and 

r +ao dk P -PHk)-iA 
*0M,fl = ^^ fc ) i- e -^( fc )-M] ^ fc ) e ^ 

is the finite temperature correction to the latter. 

It is now easy to analyze separately the three contributions (j4.15H4.l7jl . To the 
end of this section we set for simplicity m — n — 0. Then 

e B - B (0) = 6^. ( 4 - 18 ) 

which is the familiar S-B thermal energy density for a massless Hermitian scalar 
field in 1+1 space-time dimensions. 

Eq. ()4.16|) is a useful integral representation of the Casimir energy density in 
terms of the scattering matrix S(k). This representation extends the result of to 
all boundary conditions of the type (|2.5j) . In the scale invariant case for instance, 
S(k) is a constant matrix S leading to 

£ o 0M) = -A- ( 4 - 19 ) 

The corresponding finite temperature correction (|4.17|) reads 

S(x, i, j3) = ^ ^ - -S- . (4.20) 



2(3 2 sinh ( 2vr^ » bnX 
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We illustrate the case of a k- dependent S'-matrix on the example (J2.45|) of asym- 
metric boundary conditions in the vertex V. The fc-integration in (|4.16|) with the 
diagonal elements of ()2.46|) can be performed explicitly and one finds 



£ c (x, 1) 



A 



3nx 



+ 



2\ 2 e 2Xx 



3tt 



£ 



A 



+ 



A 2 e 2Ax 



Ei(-2Ax) 



Ei(-2Ax) , 



2,3 



(4.21) 



(4.22) 



127ra 6tt 

Ei being the exponential-integral function. We have plotted these functions for 
A = 1 in Fig. 2, which shows that the Casimir energy in the edge E\ behaves quite 
differently from that in the edges E 2 ^. In fact, close to the vertex along E\ one has 
an attractive force, whereas the force along is always repulsive. 
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Figure 2: Plots of £ c (x, 1) (left) and £ c (x, 2) = E c (x, 3) (right) at A = 1. 
Using the results of [13], the thermal correction (|4.17|) in this case equals 

e(x,i,p) = 

it 1 A 2X 2 e 2xX . 4 „ f(3\ 47r *\ , 

7 V " 5—2 " q o Ei -2Ax +—F[^-, e" 4 ^ , 4.23 

2/? 2 sinhf27r|J 87rx 3nx 3n P V2vr J 



e {x ,i,0 ) = -—-—W-2X x) + -F\^,e~i) , i = 2,3, (4.24) 
with 

27T(T7~ 

F(or, r) = — 2 Fi [2, 1 + a; 2 + a; r] , (4.25) 

3(1 + cr) 

where 2F1 is the hypergeometric function. Fig. 3 and Fig. 4 display the behavior of 
8(x,i,0) at fixed (3 and x respectively. 

Summarizing, we derived in this section the expectation values of some local 
observables in the Gibbs state, paying attention to the contribution of the vertex 
interaction captured by the boundary conditions ()2.5j) . In what follows we will 
address the problem of bulk interactions. 
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Figure 3: Plots of £(x, 1, 1) (left) and £(x, 2, 1) = £(x, 3, 1) (right) at A = 1. 
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Figure 4: Plots of 5(0.1, 1, (3) (left) and 5(0.1, 2, f3) = 5(0.1, 3, j3) (right) at A = 1. 



5 Non-trivial bulk interactions 



The dynamics of the classical Non-Linear Schrodinger (NLS) model on a star graph 
G is described by the equation of motion 6 

{id t + dl)^ i {t,x) = 2g\^ i {t,x)\ 2 ^ i {t,x), (5.1) 

j6l being the coupling constant. In order to avoid the presence of bound states in 
the bulk, we take below g > 0. Our main goal now will be to quantize (J5.1)) at zero 
temperature, imposing the initial conditions (|2.2I2.3|) and the boundary condition 
f)2.5p . For this purpose one must: 

(i) construct a Hilbert space Tt describing the states of the system; 

(ii) define on an appropriate dense domain D C Ti the operator valued distribu- 
tions if)i(t,x) and ij;* l (t,x), satisfying (j2.2l2.3l2.5J) and the equation of motion ()5.1j) 
with suitably defined operator product in the right hand side; 

(iii) exhibit a vacuum state Q G T>, which is cyclic with respect to the fields 

ip*' l (t,x). 

This program has been carried out for n = 1 and n = 2 in [22] and (2Sj respec- 
tively. Since the approach used there has a direct generalization to a star graph G 

6 In order to simplify the notation, in this section we fix the units in such a way that m = ^. 
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with arbitrary n, we will sketch only the main points of the construction. The first 
step is to generalize the algebra A in order to incorporate both non-trivial bulk and 
vertex interactions. Following |21]-[2E] ; we introduce the two-body bulk scattering 
matrix 7 

fli2(fc) = -r^n»®I«, (5-2) 
k + ig 

which satisfies unitarity 

R 12 {k)R 2 i{-k)=I n ®I ri , (5.3) 

Hermitian analyticity 

R* 12 (k) = R 21 (-k) , (5.4) 
and the quantum Yang-Baxter equation 

Ruih - k 2 )R 13 (h - k 3 )R 23 {k 2 - k 3 ) = R 23 {k 2 - k 3 )R 13 (h - k 3 )R 12 (h - k 2 ) . (5.5) 

In the spirit of factorized scattering [2I|-[2E| we replace ()2.19I2.20|) by 

a h (ki) a i2 (k 2 ) - Rfjl(k 2 - An) a h (k 2 ) a h (ki) = , (5.6) 
a** 1 (An) a" 2 (fc 2 ) - a*»(k 2 ) a*^(An) R^(k 2 - An) = . (5.7) 

encoding in this way the bulk interaction in the exchange relations. On the other 
hand, in order to describe the vertex interaction we introduce JTU] the so called 
boundary generators b^k), which appear in the right hand side of 

a h (An) a* 12 (k 2 ) - (k 2 ) (An - h) a n (An) = 

2n5(k l - k 2 ) 5% 1 + 27r5(k l + k 2 ) b%{k x ) , (5.8) 

being the counterpart of (j2.21|) . Due to the particular form f)5.2j) of Ri 2 , the bound- 
ary generators commute 

[#(*i),#(k)] = 0, (5.9) 
but have non-trivial exchange relations with a,i(k) and a** (A;): 

atM b? 2 (k 2 ) = R^(k 2 - An) b{\(k 2 ) R 1 ^ (An + k 2 ) a h {k x ) , (5.10) 

6£(*0 a^{k 2 ) = a* l >(k 2 ) R^ih - k 2 ) b^(h) R^(k 2 + An) . (5.11) 

The next step is to focus on the Fock representation J-'(A) of the above algebra 
A. According to (TTJ], J-'(A) is fixed (up to unitary equivalence) by the condensate 

7 We adopt the conventional tensor notation. 
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(Q , bl(k)Q), where Q is the vacuum state and (• , •) the scalar product in J-{A). In 
order to satisfy the boundary condition (|2.5|) . we require that 

(Q,b j l (k)Q) = S{(k) (5.12) 

with S(k) given by (|2.7|) . 

The final step is to define the basic structures, described in (i)-(iii) above, in 
terms of J- (A) as follows: 

• 7i and Q are the Hilbert space and the vacuum state of T(A). 

• The quantum field ipi(t,x) admits the series representation 

oo 

Ut,x) = J2(-9) m 4 m \t,x), (5.13) 



m=0 



where 



-OO 



-°° k=l 

1 = 



i E (<HX-q?t)-i £ {Pkx~p\t) 
g ;=o fe=i 

m 

II (Pfe - 9fc-i - - % - ie) 
fe=i 



(5.14) 



without summation over the upper and lower index i in the right hand side. 

• The domain T> is the finite particle subspace of J 7 (A), which is well-known to 
be dense in TC. 

It is worth mentioning that the series ()5.13|) is actually a finite sum when ipi(t, x) is 
acting on D. The coupling constant g appears explicitly in ()5.13|) and implicitly in 
a,i(k) and a* l (k), which depend on g through the bulk scattering matrix (|5.2|) . Eq. 
(J5.14j) defines a sort of non-linear Fourier transform of the type used in |2lH I3TJ] to 
solve the classical NLS model. In oder to give meaning of the equation of motion 
(|5-H) at the level of quantum fields we perform the substitution 

\^(t,x)\ 2 ^(t,x) h-> : : (t,x) , (5.15) 

(without summation over i) introducing the normal ordered product : • • • : as fol- 
lows [22] : all creation operators in such a product stand as usual to the left of all 
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annihilation operators. In view of eqs. (J5.6I5.7J) . in our case one must further specify 
the ordering of creators and annihilators themselves. We define : • • • : to preserve 
the original order of the creators. The original order of the annihilators is preserved 
if both belong to the same if) or ifi* and inverted otherwise. With this convention 
one can show [22J that 

(id t + dl){x\ , 1>i(t,x)x2) = 2 9 (Xi , ■ ^ii>* l ^i ■ (t,x)xa) , x > 0, i = 1, ...,n, 

(5.16) 

holds for any Xi,2 G D. The boundary condition (|2.5|) is also satisfied in mean value 
on T>, namely 

(Xi, [AiiP j (t,0) + Bi(d x iP j )(t,0)] X 2) = 0, Vt G R, Xi,2tV. (5.17) 

Finally, by means of (I5.6H5.11]) one can verify the commutation relations (J2.2I2.3J) . 
Since if) and if)* are unbounded operators, the subtle points in proving the above 
statements are essentially domain problems. They are faced (22 using the absence 
of poles in the analytic extension of the boundary and the bulk scattering matrices 
(J2.7J) and (J5.2J) in the upper complex fc-plane. The lack of both bulk (g > 0) and 
vertex (see eq. (J2.14Jl ) bound states is therefore crucial. 

The representation defined by (J5.13I5.14J) allows to compute all correlation func- 
tions. From the general structure of the solution one infers that the nonvanishing 
correlators involve equal number of if>i and if)* 1 . Moreover, for deriving the exact 
2n-point function one does not need all terms in the expansion (J5.13J) . but at most 
the (n — l)-th order contribution. In fact, 

(i) il (t 1 ,x 1 )r i2 (t 2 ,x 2 )) = (4 ( ;\t 1 ,x 1 )r i2{0) (t 2 ,x 2 )) , (s.is) 

(lf) h (t U X^lf)^ (t 2 , X 2 )lp* h (t 3 , X 3 )lp* U (t 4 , Xi)) = 

(4 0) (t u Xl )4^> (t 2 , x 2 )r i3{0 \t 3 , x,w im (u, x,)) + 

/(^ 0) (t 1 ,x 1 )^ 1) (i 2 ,x 2 )^ i3(1) (f3,a;3)^ (0) (^,X4)), (5.19) 

and so on. Since the vacuum expectation value of any number of {a,i(k), a* z (k), ^(fc)} 
is known explicitly, one derives in this way an integral representation for the 2n-point 
function involving n momentum integrations. The correlation functions involving 
the particle density (J4.1J) and the current (J4.2J) can be computed analogously. It will 
be interesting to extend this framework to finite temperature. 

The application of the above solution to the study of non-linear effects in the 
wave propagation on star graph quantum wires deserves further investigations. 
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6 Outlook and conclusions 



We described in the present paper a construction of quantum fields propagating 
on a star graph G with any number of edges. Our construction is based on a 
specific deformation A of the algebra of canonical commutation relations, which 
takes into account both vertex and bulk interactions. We consider the case without 
dissipation, leading to unitary theories. Two applications of the general framework 
have been discussed. The first one concerns the derivation of the expectation values 
of some observables at finite temperature. Focusing mainly on the energy density, 
we described the Casimir effect on a star graph G and derived the correction to 
the Stefan-Boltzmann law due to the interaction in the vertex V £ G. As a second 
application we solved the quantum non-linear Schrodinger equation on G. This 
example illustrates our approach at work when integrable interactions are present 
in the bulk. 

The results of the paper can be generalized in several directions. First of all one 
can increase the dimensions, setting £7, = M s and V = M s_1 with s > 1. The case s = 
2 gives rise to the so called quantum walls. In general, for s > 1 the fields propagating 
in Ei induce a non-trivial quantum field theory on V, which is no longer a point. 
The critical behavior of such kind of theories, with interactions localized exclusively 
on V, has been investigated in [3T]. On the other hand, phenomenological models 
with one edge and s = 4, aiming to construct effective theories of the fundamental 
particle interactions, have been introduced in j22] and still attract much attention. 
In this case the vertex V = M 3 is interpreted as a 3-brane confining our universe. 

Concerning the applications to condensed matter physics, it is essential to con- 
struct quantum field theory models describing the conductance properties of quan- 
tum wires. It will be interesting in this respect to extend the bosonization procedure 
and the vertex algebra construction of |3*3*1 to star graphs with any number of 
edges. The case of general graphs with N > 1 vertices, which are thus modeling 
more realistic quantum wires, represents a challenging problem as well. We are 
currently investigating some of these issues. 

Appendix 

The orthogonality relation (j2.13|) of the wave functions (j2.12|) is a consequence of 
the self-adjointness of the operator —d\ on the graph G with the boundary condition 
()2.5|) . In order check this property directly, it is convenient to introduce the T-matrix 
defined by 

S(k) =I n + iT(fc). (6.1) 
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The explicit form of T, following from (J2.61 12.71 12. 8|) . is 

T(k) = {A + \kB)- l 2\A = 2iA*{A* + ikB*)' 1 . (6.2) 

The integral over x G M + in (J2.13j) can be computed by means of the well-known 
identity 

r+ao i -I 

-i / dxe [kx = - - = ±.-md(k), (6.3) 

Jo k -\- is k 

the distribution 1/k being defined by means of the principal value prescription. One 
finds 

r+oo 

/ dxxKx; k)xi(x;p) = 5j 2<n5{k - p) + T/(k,p) , (6.4) 
Jo 

where 

v(k, P ) = -^[Ti(k)+Ti(-p)+iTmTi(-p)\ - ^ T i( k )- T i(-p)\ ■ (m 

The final step is to prove that the matrix T(k,p) vanishes. Using (J6.2j) one has in 
fact 

T(k)+T(-p) = AiA*[AA* + kpBB* + i(k — p)AB*]~ 1 A + (6.6) 
2pA*[AA* + kpBB* + i(k - p)AB*]- 1 B - 
2kB*[AA* + kpBB* + i(k - p)AB*]- 1 A , 

T(k) - T(-p) = 2pA*[AA* + kpBB* + i(k - p)AB*]- 1 B + (6.7) 
2kB*[AA* + kpBB* + i(k - p)AB*]- 1 A , 

T(k)T(-p) = -AA*[AA* + kpBB* + i(k - p)AB*]- 1 A , (6.8) 
which allow to rewrite T(k,p) in the form 

Ak 

T(k,p) = -{pA*\AA* + kpBB* + i(k - p)AB*}- 1 B - 

p 2 — k 2 

pB*[AA* + kpBB* + i(k - p)AB*}- 1 A} . (6.9) 

The vanishing of (J6.9)) now follows from the condition (J2.6|) . which implies 

[AA* + kpBB* + i(k - p)AB*]- 1 = (A* + ikB*)-\A - ipB)' 1 = 

(A* -ipB*)~\A + ikB)- 1 . (6.10) 
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